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f~^ ' Abstract 

o ' 

The linearization of complex ordinary differential equations is stud- 
ied by extending Lie's criteria for linearizability to complex functions 
/\ ' of complex variables. It is shown that the linearization of complex 

j^ ■ ordinary differential equations implies the linearizability of systems of 

partial differential equations corresponding to those complex ordinary 
differential equations. The invertible complex transformations can be 
used to obtain invertible real transformations that map a system of 
nonlinear partial differential equations into a system of linear partial 
differential equation. Explicit invariant criteria are given that provide 
procedures for writing down the solutions of the linearized equations. 
A few non-trivial examples are mentioned. 



1 Introduction 

The linearization i.e. use of point transformations to convert the nonhn- 
ear equations to hnear form of real ordinary differential equations (RODEs) 
was first introduced by Lie [7]. He provided the linearizability criteria for 
any scalar second order RODE to be point transformable to a linear RODE 
via invertible maps of both independent and dependent variables. He found 
that the RODE must be at most cubic in the first derivative and that a 
particular over-determined system of conditions for the coefficients must be 
satisfied (see e.g., [10 — 13], [15]). Tresse [2S] studied the linearization of a 
second order RODE by looking at the two relative invariants of the equiv- 
alence group of transformations, the vanishing of both of which gives the 
necessary and sufficient conditions for linearization. These are equivalent to 
the Lie conditions [15j. 

Lie's approach has been extended in many ways. The linearization of 
higher order RODEs and their systems has been discussed by various authors 
(see e.g., [IHl [191 [211 [27] ) . Further, linearization of PDEs is discussed in [5l[6l 
[13] . It is seen that for systems of RODEs, Lie linearizability criteria involves 
solving more than ten equations simultaneously that become complicated and 
messy. Also, for PDEs and their systems it is highly nontrivial to obtain the 
transformations that remove nonlinearity in the equations. In a companion 
paper, henceforth called P-I [3|, the Lie linearizability criteria is extended 
to complex functions of real variables to obtain linearizability criteria for 
systems of RODEs. 

We have extended the Lie approach of linearization of RODEs via complex 
point transformations to the complex domain and obtained corresponding Lie 
criteria for complex ordinary differential equations (CODEs). A CODE is a 
differential equation of a single complex function of one complex variable. 
The CODEs can be linearized in the same way as RODEs via complex point 
transformations. A CODE yields a system of partial differential equation 
(PDEs) after decomposing all the complex functions, variables and deriva- 
tives into their real and imaginary parts. That also yield Cauchy-Riemann 
(CR) equations which are already linear thus yield no difference for the lin- 
earization of systems of PDEs. The linearization of systems of PDEs corre- 
sponding to a CODE follows directly from the linearization of that CODE 
via real point transformations that are obtained after decomposing the com- 
plex point transformation. This is an intrinsic property of complex theory, 
which Penrose calls "complex magic" [23], that yields a non-trivial way of 



linearizing PDEs and their systems. 

The outhne of the paper is as follows. In the second section we dis- 
cuss those CODEs that are equivalent, i.e. which can be transformed into 
one another via invertible complex mappings. Lie compatibility conditions 
and some other results for linearization of systems of PDEs corresponding 
to CODEs are discussed in the third section. Examples are also discussed 
in the same section. Finally, in the fourth section we have given a sum- 
mary and discussion regarding our complex Lie symmetry (CLS) analysis for 
linearization. 

2 Equivalent CODEs 

Two CODEs are (locally) equivalent via an invertible complex transfor- 
mation if one can be transformed into the other by an invertible complex 
transformation. A CODE is a system of four PDEs, which includes the two 
linear CR equations in two unknown functions of two independent variables. 
If a CODE is equivalent to another CODE via invertible complex trans- 
formations then the system of PDEs corresponding to that CODE is also 
equivalent to the other system of PDEs. The invertible complex transfor- 
mation yields two real transformations, which can be used to transform one 
system of PDEs to another system. 
Every first-order CODE 

u'{z) = w{z,u), (1) 

can be transformed into the simplest one u' = 0, via a complex point trans- 
formation [16]. Similarly, every linear second-order CODE can always be 
transformed into its simplest form via invertible complex transformations, 
specifically into u"{z) = 0. For example, the complexified Ricatti equation 

u'{z)+u^ = 0, (2) 

is transformable to U' = by means of 

Z = z, U = {l/u)-z. (3) 

Also, the complexified simple harmonic oscillator equation 

u"{z) +u = 0, (4) 



can be transformed into U" = via an invertible complex transformation 

Z = tanz, U = u secz. (5) 

To decompose the above CODEs into corresponding systems of PDEs, write 

z = x + iy, u{z) = f{x, y) + ig{x, y), (6) 

w = wi + iw2. (7) 

The system of PDEs corresponding to a first-order CODE is 

fx + 9y = wi{x,yj,g), g^ - fy = W2{x,yJ,g), (8) 

which can be transformed into 

F^ + Gy = 0, G^-Fy = 0, (9) 

by an invertible real transformation derived from the complex point transfor- 
mation. Similarly, the system of PDEs corresponding to a linear second-order 
CODE is 

fxx- fyy + '^Qxy = Wi{x, yj, g, h, I), (10) 

gxx - gyy - '^fxy = W2{x, y, f, g, h, I), (11) 

where both wi and W2 are two such real functions that do not give rise to a 
nonlinear system. The above system can be transformed into 

Fxx - Fyy + 2G'xy = 0, (12) 

Gxx — Gyy — '^Fxy = 0. (13) 

via an invertible real transformation that is obtained from a complex trans- 
formation. 

Examples 

1. The complexified Ricatti equation is equivalent to 

fx + gy = -f^ + g\ (14) 

g.-fy = -2fg. (15) 



If we apply the following transformation 

X = x,Y = y, (16) 

X, G= —- — --y, (17) 



we get 

Fx + Gy = 0, (18) 

Gx-Fy = 0. (19) 

Notice that a simple complex transformation for the complexified Ricatti 
equation yields a non-trivial real transformation (16) and (17) that trans- 
forms the system of PDEs (14) and (15) into its simple analogue (18) and 
(19). This remarkable feature of complex theory is of great significance for 
the linearization of PDEs and their systems. 

2. The system 

Jxx ~ Jyy + ■^9xy = ~/; l^Uj 

dxx ^ 9yy ^ ■^Jxy = ~fi'; l^-'-j 

can be transformed into 



XX — ^YY 



2GxY = 0, (22) 



Gxx — Gyy — '2FxY — 0, (23) 

via the real transformations 

_ (1/2) sin 2x _ (l/2)sinh2t/ 

cos^ X + smh y cos^ x + smh y 

/ cos X cosh y — (7 sin x sinh y „ f sin x sin y + (? cos a; cosh y 

^ ^ 2 ^ ^T2 ' ^ = 2 ^ ^T2 • (2^) 

cos^ X + smh y cos^ a; + smh y 

Thus, the system of PDEs corresponding to a CODE can easily be trans- 
formed into its simplified (linear) form by invoking complex transformations. 
It would have been very difficult to guess or calculate the real transforma- 
tions (24) and (25) that transform system (20) and (21) into system (22) and 
(23) without this formalism. It is an extraordinary characteristic of complex 
magic. 



Lie conditions in the complex domain can be derived by replacing the 
real variable x in r-CODE in P-I with a complex variable z. A scalar second 
order CODE which is at most cubic in its first derivative 



"/'~^ _ /1/'~ „,A„ ,'3 I D/'~ „A„ ,'2 



u (z 



A{z, u)u^ + B{z, u)u'^ + C{z, u)u' + D{z, u), (26) 



where A, B, C and D are complex valued functions, is linearizable according 
to Lie theorem. We do not re-state Theorem 1 of P-I as all the conditions (1- 
9) remain the same with only the change that the real variable x is replaced 
by a complex variable z. In the subsequent section we will decompose (26) 
into system of PDEs by using (6) and by taking h and / as real and imaginary 
parts of u' . 

3 Lie Conditions for Systems of PDEs 

According to Lie's theorem the necessary condition for a second order 
CODE to be linearizable is that it is at most cubic in its first order derivative 
implies the necessary conditions for system of PDEs corresponding to that 
CODE i.e. the system must be at most cubic in first derivatives together with 
certain constraints on coefficients. The real transformations for linearization 
of a system of nonlinear PDEs can be obtained by decomposing complex 
transformations that linearize a CODE. The general form of a system of 
PDEs corresponding to (26) is given by 

(27) 
g.c.c-gyy-2f^^y = A\3h'^l-l^)+A\h^-3hl^)+2B^hl+B^{h^-l^)+C^h+C^L+D^. 

(28) 
where all the coefficients A\ B^, C* and D^ are functions of x, y, f, g. 

Now we derive Lie conditions for systems of PDEs corresponding to 
a second order CODE by decomposing results in Theorem 1 in P-I into real 
variables by invoking (37) and (38) and by taking p as a complex function of z. 

Theorem 1. The following statements are equivalent. 

1. The system of PDEs (27) and (28) is linearizable via real transformations; 



2. The coefficients in (27) and (28) satisfy 

3Ai, - ?,Aly + QAly + ?>C^Al + ?>C^Al - ?>AlC^ + SC^Aj - SA^D^ - 3^^^^ 
^D^A) - SD^Al + SA^Cl + 3A^C^ - SA'C^ + 3A^C^ + C]^ - 0]^ + 20}^- 
QA^D) - QA^D] + QA^D) - QA^D] + B^C] + B^Cl - B'^CJ + B^C^ 



2B'Bl - 2B'Bl + 2B'Bl - 2B'BI - 2Blf - 25j^ - 2Bl^ + 25^^ = 0, 



(29) 



d-/i^^ 0-/l^^ b-/l^^ ~r oL/ -/l^ ~r oC- ./l^ ~r o/l^C' oC- ./l^ di.^' ./Ir oLJ -/i^" 

SDMJ + 3DMJ + SA^Cl + 3A^C^ + 3A^C^ - SA^Cl + C}^ - C]^ - 2CJg- 

GA^D} - QA^D] - 6A^DJ + GA^DJ + B'^CJ + B^C'^g + B^CJ - B^C]- 

2B'Bl - 2B^Bl - 2B^Bl + 2B^Bl - 2B% + 25^^ + 25^^ - 25^^ = 

(30) 



QD^Al + 6DMJ - QD'^Al + GD^Aj - 3D^5} - 3D^52 + 3D25J - 3D^Bl + 

SA'Dl + 3AiZ}J - SA'Dl + SA'Dl + B^ - B^ + 2B% - 2Clj - 2Cj^- 
2C% + 2Clg - 3^1^) - W^D] + 352Z}2 _ 352^1 ^ 3^1^ _ 3^1^ ^ g^2^^ 

2^1^) + 2C^CJ - 2C^C} + 2^2^] - C^E^ - C^Bl + ^^5^ - C^B] = 0, 

(31) 

6D=^Ai + QD^Al + QD^Al - QD^A\ - 3D^BJ - 3D^B^g - 3D^BJ + 3D^Bl + 

3A'Dl + 3A'Dl + 3A'Dl - 3A'dI + B^ - B^ - 25^^ - 2C',f + 20]^+ 

2Clg + 2Clg - 3B^D) - 3B^Dl - 3B^DJ + 3B^Dl + 3DJf - 30^^ - QD}g+ 

2C^CJ - 2C^Cl + 2C^C} - 20^0] - C^Bl - C^B^ - C^B^ + C^B] = 0; 

(32) 

3. The system of PDEs (27) and (28) has four real symmetries Xi, Yi, X2 
and Y2 with 

Xi = P1X2 - P2Y2, Yi = P1Y2 + P2X2, (33) 

for nonconstant pi and p2 and they satisfy 

[Xi, X2] - [Yi, Y2] = 0, [Xi, Y2] + [Yi, X2] = 0, (34) 



such that a point transformation {x,y,f,g) — > {X,Y,F,G), which brings 
Xi, Yi, X2 and Y2 to their canonical form 

d d d d d d 

reduces the system to the hnear form 

Fxx - Fyy + 2GxY = W,{X,Y), (36) 

Gxx - Gyy - 2FxY = W2{X, Y); (37) 

4. The system of PDEs (27) and (28) has four real symmetries Xi, Yi, X2 
and Y2 with 

Xi = P1X2 - P2Y2, Yi = P1Y2 + P2X2, (38) 

for nonconstant pi and p2 and they satisfy either 

[Xi, X2] - [Yi, Y2] ^ or [Xi, Y2] + [Yi, X2] ^ (39) 

such that a point transformation {x,y,f,g) — > {X,Y,F,G), which brings 
Xi, Yi, X2 and Y2 to their canonical form 

reduces the system to the linear form 

Fxx - Fyy + 2Gxy = HW,{X, Y) - LW^iX, Y), (41) 

Gxx - Gyy - 2Fxy = HW2{X, Y) + LWi{X, Y). (42) 

where 

H = Fx + Gy, L = Gx-Fy; (43) 

5. The system of PDEs (27) and (28) has four real symmetries Xi, Yi, X2 
and Y2 with 

Xi ^ P1X2 - P2Y2, Yi ^ P1Y2 + P2X2, (44) 

for nonconstant pi and p2 and they satisfy 

[Xi, X2] - [Yi, Y2] = 0, [Xi, Y2] + [Yi, X2] = 0, (45) 



such that a point transformation {x,y,f,g) — > {X,Y,F,G), which brings 
Xi, Yi, X2 and Y2 to their canonical form 

Xi = ^' Yi = — , ^2 = ^,Y, = — (46) 

reduces the system to the system of PDEs which is at most cubic in all its 
first derivatives; 

6. The system of PDEs (27) and (28) has four real symmetries Xi, Yi, X2 
and Y2 with 

Xi ^ P1X2 - P2Y2, Yi ^ P1Y2 + P2X2, (47) 

for nonconstant pi and p2 and they satisfy either 

[Xi, X2] - [Yi, Y2] ^ or [Xi, Y2] + [Yi, X2] ^ (48) 

such that a point transformation {x,y,f,g) — > {X,Y,F,G), which brings 
Xi, Yi, X2 and Y2 to their canonical form 

d d d d d d 

d d d d 

reduces the system to the linear form 

X{Fxx - Fyy + 2G'xy) - Y{Gxx - Gyy - 2Fxy) = ai{H' - 3HL')- 
a2{3H^L - L") + 61 (/J2 _ L") - 2b2HL + -—^-^-{3{al + 4) + (6^ - bl)ai+ 

o( (2-1 ~r 0^0 ) 

2bib2a2}H - {2&i&2ai - a2(&? - bl)}L + ^ ^ A ^iai + 6202) + 

{{bl - 3bibl){al - a\) + 20102(36162 - 63), (51) 



27(a2 + al) 

Y{Fxx - Fyy + 2Gxy) + X{Gxx - Gyy - 2Fxy) = a2{H^ - 3HL^) + 
a,i3H^L - L") + b2{H^ - L^) + 2b,HL + / {3(a^ + aj) + {bj - 6^)01+ 

2616202}^ + ^ 2^(2^1^201 - 02(61 - 62)}i^ + ^, 2 , 27(^2^1 - 6102) + 

3[ai + 02) 3(ai + 02) 

^ {(36262 -63)(a2-a2)- 2(6? -36i62)aia2}. (52) 



27(af + al 



The invertible real transformations 

^ = F + — ^ j^iihai + b2a2)X - (6201 - ha2)Y}, (53) 

3(af + 02) 

1 
Y = G+ ^i 2 , 2x {(^iQi + ^2a2)>' + (&2ai - 6102)^}, (54) 



2^^ ^^"^3(a2 + a2) 



F = ^(F^ - G^) + 77772^T2v[{(^i«i + ^202)^ - (6201 - hia2)Y}F- 

1 



{(6iai + 6202)^ + (6201 - hia2)X}G] + t^7^^^t^[{(&? - ^2)^ - ^2) + 
46i62aia2}(X2 - F^) _ 2XY{{2hMai - 4) - 2aMbl - bl)}] + 

^ -{a^{X^-Y^)+2a2XY}, (55) 



2(af + al 
G = FG+ —^ 2^[{ibiai + b2a2)X - (6201 - bia2)Y}G+ 

{{ha, + h2a2)Y + (6201 - ha2)X}F] + ^\ [2XY{{hl - hl){al - al) + 

461620102} - (X2 - Y^){2hb2{al - al) - 2aia2{bl - bl)}] + 
1 



-{2XFai-a2(X^-F^)}. (56) 

^1^1 + a2J 

transform (51) and(52) into the system of hnear PDEs 

Fxx - Fyy + 2G^y = 0, (57) 

'^XX ~ ^YY ~ '^FxY = 0. (58) 



Examples. Now we discuss some illustrative examples. 
1. Consider a nonlinear CODE of the form 

u" + 3uu' + M^ = 0, (59) 

which is linearizable as it satisfies the Lie complex conditions. The set of two 
non-commuting CLSs is 

Zi = |-, Z2 = z^-u^. (60) 

Bz Bz Bu 

10 



Note that Zi 7^ p{z,u)Z2. We invoke condition six of Theorem 1 to find a 
hnearization transformation. The complex point transformation that reduces 
the symmetries (60) to their canonical form is 

Z = -, U = z+-, (61) 

u u 

and (59) reduces to 

Zf/" = -f/'3 + 6f/'2-llf/' + 6, (62) 

by means of the transformation (61). Equation (59) linearizes to U" = via 
the complex transformations (25) by placing a = —1, 6 = 6 in equation (34) 

in P-I. That is 

1 ~ z'^ z 
Z = z--, U = — --. (63) 

It may be seen that the above transformation is a proper complex trans- 
formation thus one does not get a contradiction as was in P-I. Further, the 
procedure that seems to be analogous to analytic continuation do not arise 
here. The system of PDEs corresponding to (59) is 

f.. - fyy + 2g,y = -3{fh - gl) - {f - 3//), (64) 

9xx - Qyy - '2fxy = Sigh + fl) - (sfg - /), (65) 

which admits the RLSs 

Xi = |-,Yi = |-, (66) 

ox ay 

Y d d d d 

ox oy oj og 

^^ d d d ^ d , . 

Y2 = i/T- - X— - (?— + /— . 68 

ox oy OJ og 

These generators satisfy (47) and (49). The transformation 

fl ^gl p^ gl 

P = ^ + -J—2^^ = y--rr—2^ (^°) 

J +9 J +9 



11 



transforms (64) and (65) into 

X{Fxx — Fyy + SGxy) — Y{Gxx — Gyy — '^Fxy) = —H + ?>HL — 

6iJ2-lliJ + 6, (71) 

X{Gxx - Gyy - 2Fxy) + Y{Fxx - Fyy + 2Gxy) = L^ - SH^L+ 

12HL - IIL, (72) 

where H and L are given by (43). The above system can further be reduced 
into very simphfied hnear PDEs by using 

^ = ^^"72^ — 2^yf-^9)- (75) 

The hnear PDEs are 

hx - Fyy + '^GxY = 0, (76) 

Gjix - Gyy - 2-^xy = 0. (77) 



2. Consider a second order nonhnear CODE with an arbitrary function 

w{z) 

uu" = u'^ + w{z)u^. (78) 

The above CODE admits two CLSs of the form 

Zi = zu—, Z2 = u—, (79) 

and thus 

[Zi, Z2] = and Z2 = -Zi. (80) 

By using the complex transformation 

Z = - and [/= -logM (81) 

z z 

12 



(78) can be reduced into a linear CODE 

U" = ^,w{^). (82) 

The system of PDEs corresponding to (78) are 

fifxx-fyy+2gxy)-9i9xx-9yy~2f^y)=h'^-l'^+wi{f-g'^)~2fgw2, (83) 

fi9xx-gyy-'2fxy) + ifxx-fyy+2gxy)9=2hl+2Wifg+W2if-g'^). (84) 

Invoking the transformation 

X = "^ Y = ~^ fSS) 

F = -^/-xHf + g')+yt^n-\^)], (86) 

X -\- y z J 

G = ^^[xtan-i(f ) - I Hf + /). (87) 

x^ + y^ J 2 

Equations (83) and (84) reduce to the hnear system of PDEs 

Fxx-Fyy+2Gxy= ^, ] ^^ [{X'-3XY')wi-{Y^-3X^Y)w2], (88) 

X"^ + Y "^ 

Gxx-GYY-2FxY=—^[{X^-3XY')w2+{y'-^X'y)wi], (89) 



where 



V Y X —Y 

^1 = ^i( j^2 + y2 ' x^ + y^ )' ^2 = ^ x^ + Y^ ' x^ + y^ ^- ^^°^ 



3. Consider the nonhnear CODE 

m" = -(m' + m'=^). (91) 



This equation has CLSs 



Z, = iA, z, = -^, (92) 

z (72; 2; (72; 
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The complex transformation 

U = -z^ and Z = u, (93) 

transforms (91) into a linear CODE 

[/" + 1 = 0. (94) 

The system of PDEs corresponding to (91) is 

/.. - fyy + '^g.y = ^^^[xh{l + h^ - Sf) + yl{l -P + 3h% (95) 

1 

x^ + y 
and admits CLSs 



9.. - 9yy - 2Uy = -z^—^lxKl - /' + 3h^) + yh{l + h^ - 31% (96) 



X _ ^ ^ y ^ Y - ~y ^ X d 

^2 _|_ y2 Qrj. 2;2 jf- y2 Qy^ x^ + y"^ dx x"^ + y"^ dy^ 

X _ f^ ^ fy d yg d gx d 

J,2 _|_ y2 g^ ^2 _|_ y2 Qy ^2 _|_ ^2 g^ rj,2 _|_ ^2 ^^ ' 

Y _ ^a: (9 ^y ^ /y "g jx d 

rj.2 _|_ ^2 ^2; x^ + y^ 9t/ a;^ + y"^ dx x^ + y'^ dy 

The above RLSs correspond to the condition five of Theorem 1. The real 
transformation that linearize system (95) and (96) is 

X = f,Y = g, F=]^{x^-y^),G = xy. (100) 

The linearized system of PDEs is 

Fxx - Fyy + 2Gxy + 1 = 0, (101) 

Gxx - Gyy - 2FxY = 0. (102) 



4. Consider the nonlinear CODE 

u" = 1 + (u' - zfw{2u - z^) (103) 

14 



which admits the CLSs 



d d „ (9 2 "^ 



oz ou oz ou 

Equation (103) reduces to a hnear CODE by using the complex transforma- 
tion 

Z = 2u-z^ andU = z. (105) 

It becomes 

U" = ^U'w{Z). (106) 

The system of PDEs corresponding to (103) is 

fxcc - fyy + ^g^y = I + {{h - xf - (/ - yf}w, - 2{h - x){l - y)w2, (107) 
9xx - 9yy - '2fxy = 2(/i - x){l - y)wi + {{h - xf -{I- yf}w2, (108) 

which admits RLSs 

^r 9 d d ^^ d d d ,^^^, 

ox of Og og Oj Oy 

which satisfy condition six of Theorem 1. The transformation 

X = 2f-x^ + y^,Y = 2g- 2xy, (112) 

F = x, G = y. (113) 

yields system of linear PDEs 

Fxx-Fyy + 2Gxy = --{Hwi-Lw2), (114) 

Gxx - Gyy - 2FxY = -^iHw2 + Lw^), (115) 

where 

w{z)=wi+iw2, H = Fx + Gy, L = Gx-Fy. (116) 
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4 Conclusion 

Finding exact solutions of nonlinear DEs is difficult. Linearization of 
DEs is a method to convert the original nonlinear DE to an equivalent linear 
DE so that exact solutions can be constructed directly. But it requires the 
existence and construction of such transformations (point, tangent, contact or 
complex) that transforms nonlinear DEs into linear DEs. For scalar ODEs Lie 
provided a nontrivial construction of point transformations via symmetries, 
but it is difficult to extend the Lie conditions for systems of ODEs and for 
PDEs (see e.g., [SI El [13], [15 — 19]). Finding transformations that map 
systems of nonlinear ODEs and PDEs into linear systems is not only tedious 
but also highly nontrivial. 

We have applied CLS analysis to find non-trivial ways of reducing order, 
linearizing and solving certain systems of ODEs and PDEs. Several known 
results of classical symmetry method were studied in the complex domain. 
Various results for systems of PDEs were obtained, e.g. the use of CLS 
analysis to study CLSs for CODEs and respective real Lie symmetries for 
systems of PDEs corresponding to CODEs, discussed in [1]. Variational 
problems involving complex Lagrangians of CODEs results in Lagrangians 
of systems of PDEs corresponding to those CODEs and double reduction 
in order of CODEs via complex Noether symmetries, which are studied in 
[2]. Further, variational problems for r-CODEs and their use in constructing 
Lagrangians for systems of RODEs are discussed in [1] . 

We have looked at linearizability criteria for second order CODEs and 
obtained analogous Lie conditions for CODEs and applied them to the lin- 
earization of systems of nonlinear PDEs in this work. The corresponding 
statement for a system of two PDEs (in fact four i.e. by incorporating linear 
CR-equations) associated to a second order CODE was presented. Examples 
were given for the linearization of systems of PDEs including the construction 
of the point linearizing transformations. 

In this paper the complex function is a function on the complex plane 
that requires CR-equations to hold for complex differentiability. Thus, we 
mapped complex solutions of CODEs into complex solutions via complex 
transformations. Since, CR-equations are linear this makes no difference in 
the linearization of systems of PDEs corresponding to CODEs except that 
all the symmetries are analytic i.e. their coefficients satisfy CR-equations. 
These form the conformal subalgebra admitted by systems of PDEs. It is 
hoped that further classification of systems of PDEs in [1] with respect to 
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conformal algebras admitted by them can be achieved. 

We have seen that hnearization theorems can be carried over from the 
real case to the complex case with remarkable results for the linearization of 
systems of PDEs. Thus we were able to linearize those systems of nonlinear 
PDEs that correspond to some CODEs. There are several possible ways of 
using CLS analysis. It would be of great interest to extend the geometric 
proof of linearization for RODEs [T0l|2Q] to CODEs. The analysis can also be 
used for classification of those systems of nonlinear PDEs that correspond to 
CODEs. An extension of the Lie linearization results for systems of quadrat- 
ically and cubically semi-linear CODEs in [191 EQ] , in the complex domain 
may yield useful and highly nontrivial results for systems of PDEs. Further, 
results of linearization for third and fourth order CODEs may be obtained 
by using results for third and fourth order RODEs in [TH [211 122] • 

It is expected that the Lie criteria for linearization of two CODEs can be 
found in a similar way to the above, which can then be used to linearize four 
nonlinear PDEs corresponding to a system of two CODEs. A system of two 
second order RODEs admits 5, 6, 7, 8 or 15 RLSs and the maximal symmetry 
algebra is s/(4, R ) for the simplest system [27j. Thus, it may be argued that 
the system of two second order CODEs also admit 5, 6, 7, 8 or 15 CLSs. These 
result in 10, 12, 14, 16 or 30 RLSs which form the subalgebra of symmetries 
admitted by four second order PDEs that correspond to two second order 
CODEs. Further, maximal symmetry algebra s/(4, (D ) is admitted by the 
simplest systems of two second order CODEs. 

It is also hoped that Penrose's "complex magic" can be further extended 
to ^^hypercomplex magic " by using hypercomplex variables (e.g. quaternions, 
octonions, Clifford or Grassman variables) that should provide hypercomplex 
Lie symmetries. These can then be used to reduce the order, linearize and 
solve a large class of systems of four PDEs by using a scalar ODE in hyper- 
complex variables. It is further expected that the non-commutative behavior 
of hypercomplex variables will yield a non-trivial generalization of Lie's work. 
Further, use of the hypercomplex Lie symmetries in geometric calculus ^ 
may result in non-trivial physical implications. 

In P-I, we used CLS analysis to construct linearizability criteria for cer- 
tain systems of nonlinear RODEs by considering r-CODEs. An r-CODE is a 
complex ODE in a single complex function of one real variable that yields a 
system of two RODEs. The linearization of r-CODEs directly gives lineariz- 
ability criteria for systems of RODEs corresponding to those r-CODEs [3]. 
Again, the use of hypercomplex variables should yield linearizability criteria 
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for systems of four ODEs. The possibilities to extend Lie's programme of 
linearizing nonlinear differential equations, and determining which equations 
can be so reduced, seem endless. We have only touched the tip of the iceberg. 
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